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t*"** ! Instanton liquid in heated and strongly interacting matter is studied using the variational princi- 

ple. The dependence of the instanton liquid density (gluon condensate) on the temperature and 

fNj ■ the quark chemical potential is determined under the assumption that, at finite temperatures, 

the dominant contribution is given by an ensemble of calorons. The respective one- loop effective 
quark Lagrangian is used. 

<n : 

In current studies of strong-interacting matter under extreme conditions, primary attention is 
; focused on a description of its phase state at given temperature and chemical potential. For def- 
i ■ initeness, we consider that T is the temperature of quarks and \i is the quark chemical potential 
(it is assumed that gluons are in thermodynamical equilibrium with quarks). However, there is no 
yS^, approach making it possible to describe main features of the expected phase diagram of quark-gluon 
\ matter at least qualitatively. 

In the present study, we argue that the instanton liquid model of the QCD vacuum [I] can shed 
' light on some important features of a full picture. It is frequently noted that this model offers a useful 
tool for obtaining phenomeno logically plausible estimates in spite of the fact that it is poorly justified 
q ! because the typical size of an instanton is not properly fixed. As of now, this fact is considered as 
^ | inessential because a connection has been revealed between limitations on the instanton size due 
O t° repulsion [2] and generation of mass of the gluon field in the framework of the quasi-classical 
approximation [3]. The latter mechanism is a more general property of stochastic gluon fields than 
the former one. We will discuss this question later. Here we assume that the problem of instanton 
• i—i ! size is solved in one of the following scenarios: self-stabilization of the saturating ensemble [2],[1], 
freezing of the coupling constant [3], or influence of the confining component [6]. In the present 
c3 ; study, primary attention is focused on a plausible qualitative model describing a behavior of the 
gluon condensate. 

In the beginning, we recollect the variational principle proposed in [2] and the method of determi- 
nation of the size of pseudoparticles and the density of the instanton liquid and introduce notation for 
further considerations. In the model of instanton liquid describing the QCD vacuum, it is assumed 
that the leading contribution to the QCD generating functional is given by the background fields 
representing superposition of instantons in the singular gauge: 

A K x ',l) = - ^ ab Vb^ a v {y) , a„{y) = — — 5 -| , y = x-z, fi, v = 1, 2, 3,4 . (1) 

g y +p y 

where p is the size, u is the matrix of color rotation, and z is the position of the center of a 
pseudoparticle (in the case of anti-instanton, the 't Hooft symbol should be replaced as follows: 
?7 — ^ 77) . This being so, the QCD generating functional takes the form 

00 1 N . 00 1 JV . 

y = e 4 n ** « ^ Uint{i) = e 4 n ** e " £(7) > ( 2 ) 
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E{i) = (3 U mt (l) ~ E ^d( Pl 
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d{p) = 4 P 2Nc e~ m , (3) 



P 

is the instanton size distribution [7]; d^i = dzi duji dpi, and 

9 

is the action of a single instanton, where (A = = 0.92Ap y ) CV c , depends on the renormalization 

scheme and, in the case under consideration, is given by CV ~ 4-66 exp( 1 .68A^) ^ _ 

& 7r (AT C — 1)!(AT C — 2)! 

- - ^ c ^ - —£ . We assume that Nj=2 here because the leading contribution to renormalization 
comes from hard massless gluons and quarks. The auxiliary function 

p = -b ln(Ap) , 

is evaluated at the scale p defined by an average size of pseudoparticles, Umtij) is considered assuming 
pair interaction dominance. Its contribution has the form [2] 

J dw x dto 2 dz x dz 2 £^(71,72) = V £ 2 p\ p\ , 

where £ 2 = 27 ' ^ N ^ ^ . The factor (3 that appears in the exponent in formula ^ is also evaluated 

at the scale of an average size of pseudoparticles p. Assuming that the instanton liquid is topologically 
neutral, we do not introduce notation to distinguish between instantons and anti-instantons, A" 
denotes the ovarall number of pseudoparticles in volume V. 

Since the interaction is independent of coordinates or orientation in color space, it is natural to 
calculate the generating functional Y on the basis of the effective one-particle distribution function 
p(p), which can be determined from the solution of the variational problem 



00 yN N . 00 yN N , 

J dp, P(P) = E W n Jdli e~ E ^ , (4) 

JV=1 i=l 7V=1 i=l 

^1(7) = ~E ln P(pi 



where the factor V N in (j3J) is isolated in order that the result be expressed in terms of the respective 
density and convenience in interpretation of the function p(p). With regard to convexity of the 
exponential function, the generating functional (j2J) for every fixed A" partial contribution can be 
estimated using the approximating inequality 

Y'>Y a = Y{ expi-iE-Et)) , (5) 

where an average over approximate ensemble is implied. In the case under consideration, the average 
of difference (E — E±) is given by: 

1 1 r N 

(E-E 1 ) = — —J J] dn \p U mt -J2^d( Pl ) + J2^P(Pi)) e E ' nMft) = 

= ^ (/ rfp ^ ln djp} + 2p a V I dpldp2 ^ p i p 2 P(pi)p(p2)j , 
where po = J dp p(p). 
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Variation of the functional (E — Ex) with respect to p(p) results formally in the equation p(p) = 
e~ 4 d(p) e~ n ^ 2p2p2 (where n = N/V is the density of the instanton liquid). Here an unwanted factor 
of e _1 , emerges. It can be excluded due to the fact that the approximate functional Y a is independent 
of the constant factor of C that can be added to the expressionfor p(p). For convenience, we set 
C = e, and therefore, arrive at 

p(p) = d(p) e-'^W . (6) 
Substituting this solution to the approximate functional, we obtain 



V N ^ 
N\ 



■ 
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(P 2 ) 



Defining suitable parameter v the integral for determination p can be represented in the form 



p = A 4 J dpA C Nc f"° (pA) 
From the comparison of which with formula ([6]) we obtain 

K 2 n7 ■ 
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(7) 
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Provided that v is known, this formula offers a relation between the average instanton size and the 
density of the instanton liquid. To find this relation, we consider the equation 



J dp p 
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J dp p 
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p 2 v~ x r(^ + r 
r(^) 



It should be noted that the factor 



which gives v — and therefore, p = A 4 Cn c P 2Nc ^jj) 

of two in the denominator of this expression stems from the integration measure 2pdp, which, in its 
turn, emerges in transformation to the Gaussian integral with respect to p squired. This factor was 
omited in [2J; however, this fact has no noticeable consequences. The reason is that the parameter 
A is determined from a fit to some observable, for example, to the pion decay constant. In so doing, 
everything is governed by a choice of scale. Moreover, it should be remembered that the instanton 
liquid model is merely a rough approximation. From the above, we derive an approximate expression 
for the functional as follows: 



Y a = exp <^ —N 



[ln(n/A 4 ) - 1] + iVm 



c N 3 2N m 2 v)- u/2 



(9) 



Now we find the value of n at which the argument of the exponential approaches its maximum. To 
do this, we should solve the equation 



^ + ln(n/A 4 )+ln 
From the relation ([8]) we obtain 



r(u) 



+ n 



2N c d(3 
/3 dn 
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I d3 ldn 4 
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On the other hand, -Jj= = — ~, -^L = We represent the derivative of P with respect to the density 



in the form -M- 
dn 



dp idn 
lip' Hp 1 



and obtain 
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Thus we derive the expressionfor the instanton liquid density 



n/A 4 



r(u) 



2 
i/+2 



cxp 



' 4iV c (3 
u + 2/3 



v 



u + 2 



v 



2/2/3 



v 



(12) 



The contribution of the derivatives of the functions (3 and (3 with respect to the density was disre- 
garded in [2]. This contribution compensates for the above-mentioned factor of 2 though, as was 
noted above, this is not essential. The obtained formula for the instanton liquid density by itself 
does not provide a solution to the problem because it remains to solve the transcendental equation 
([H]) in p, where the function (3 involves the logarithm of p. To solve this equation, it is convenient to 
reformulate the problem without resort to the explicit formula (fT2l for the instanton liquid density. 
By definition of the function /3, the action of an isolated pseudoparticle must be positive. This gives a 
limitation to the maximum size of an (anti-)instanton as follows: pAC^ < 1 (actually, pA < 1). Now 



we can solve the transcendental equation ( flOl) . by bisection of the segment. In so doing, a stationary 
value of p is determined at each step and the respective instanton liquid density is determined from 
equation ([H]). In the calculation of the generating functional, the contributions of the type {pA) 2v are 
used rather than the expression for the instanton liquid density. 

Now we modify the variational principle in order to extend our description to the case of finite 
temperatures. For this purpose, we employ calorons — solutions of the Yang-Mills equations periodic 
in the Euclidean time. The background field should be replaced by a superposition of calorons and 
anti-calorons as follows [8]: 



n 



uJ ab Vb„u ^hill, 



1 + 



np 2 T 



sinh(27rrT) 



(13) 



r cosh(27rrT) — cos(27ttT) 



where T _1 is the period of the caloron, r = \x — z\ is the distance from the center of the caloron z in 
three dimensional space, and r = x 4 — - is the respective interval of "time". As the temperature 
tends to zero, such solutions go over to (anti-)instantons in the singular gauge. Yet another modifi- 
cation of the variational principle is the replacement of the distribution ([3]) in the instanton size by 
the function 

1 



d(p,T) 



/3^exp[-/3(p) 



A Nc T*p*} 



(14) 



r 2 accounts for the additional contribution to the action 



where the coefficient = ^ (^-N c — l) 
of each individual pseudoparticle. It provides an approximation to a more exact expression 



d(p, T) = d(p, 0) exp < — 



2 2 (N c + N f /2) 47TV 
9 3 g* 



+ 12 A(npT) [1 + (N c - N f )/6) 



, (15) 



constructed from the respective determinants [9]. For our purposes it is sufficient to say that the 
function A(npT) is determined by a shape of the pseudoparticle ([TBI . This function was studied in 
the cited work; however, we do not use it in the present article. It should be mentioned that the 
expansion up to the terms of the order T 2 can be used as an approximate expression for the function 
A(jipT) because, within the accuracy of the variational principle, only the terms up to order p 2 
should be kept in the argument of the exponential in formula f Tl4l) . The first term in formula ( fl5l) 
is represented as a product of two factors; each factor was interpreted in [9]. The first factor is 
the square of the electric mass, that is, the temporal component of the gluon polarization tensor 
evaluated at the zero energy and momentum. It has the form 



n 44 (w = o,p = o) 



g 2 T 2 



(N c + N f /2) 



(16) 
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The remaining components being equal to zero at zero energy-momentum. Therefore, the magnetic 
mass vanishes. Note that the one-loop quark and gluon contributions to the polarization tensor 
are taken into account [10J, the resulting sum being rearranged in order that the quark and gluon 
contributions in the medium sum up to a finite value. This, formally, gives rise to a generation of 
the mass of the gluon field. The second factor is the integral of the square of the fourth component 
A± of the field in formula (fT3l) 

4ir 2 p 2 



dy At(y)At(y) 



9" 



(17) 



It is independent of the temperature [TT]. It is seen that one can take into account only one- loop 
contribution ^m 2 el A\A\ to the Lagrangian of the gluon field and neglect other corrections. It was 
demonstrated [3] that the term U int describing the interaction of pseudoparticles can be brought in the 

Thus the interaction term also describes generation of 



form T^m 2 A^A 1 ^, where m 2 



9ir 2 n p 2 



N 2 -l 



the mass of the gluon field in the instanton-anti-instanton medium in quasi-classical approximation. 
This being so, chromoelectric and chromomagnetic fields are screened equally well provided that 
the instanton liquid density is not equal to zero. It was shown that screening is a consequence of 
stochastic character of the ensemble of gluon fields being unrelated to a specific instanton solution 
of the type ([T]) or details of the repulsion mechanism responsible for stabilization of the ensemble 
[2]. An application of these considerations to the (anti-)instanton solution (TjQ) leads precisely to the 
formula for U^. It turns out that, in the caloron ensemble, screening of chromomagnetic fields and 
the interaction term depends only weakly on the temperature. However, the anisotropy is negligible 
small and the interaction term coincides with that obtained for the (anti-) instanton solution. First 
it was found in [TT], where the instanton liquid was studied at non zero temperature. 

The one- loop contribution of Plank gluons is proportional to N c (see formula ( TT6|) ) and does not 
vary as the chemical potential becomes different from zero. On the other hand, it is known that 
the one-loop fermion contribution in the medium can be calculated exactly. It has no dangerous 
singularities [12], [13]. The "temporal" component of the polarization tensor generated by a quark 
of definite flavor has the form 



n{ 4 (fc 4 ,u;) 



k 2 



9' 



7T 2 U 2 



pu 



dp p 2 



n.,. 



As 2 
1+ v 



8pk 



k 2 (k 2 
In 



2puj) 2 + 4e 2 p k 2 



arctan 



8pu Epki 



Aefti - Ap 2 u 2 + k 4 



(k 2 - 2pu) 2 + Ae 2 k 



where u 



n, 



\k\ 



, k 2 - 



n 



+ Bp 

= (e t + 1 



(m + p 



2U/2 



where m - is the quark mass, 



n„ 



n„ + nj. 



tensor takes the form 



2 r c 
U f (h,u) = g 2 — \ 

TT JO 



dp p 2 



After summation over all components, the polarization 



n.„ 



1 + 



2m 2 - k 2 



8pk 



In 



(k 2 + 2puf + 4fflj 
{k 2 - 2pu) 2 + Ae 2 p k 2 



18) 



It is seen that, at k^ = 0, and small values of u, the first term (that is, unit) gives the dominant 
contribution to the gluon mass. The spatial components are negligibly small. In particular, at u = 
we obtain 

l44 (0,0) = s 2 ^ (19) 



^(0,0) 



TT" 







and at T = we arrive at 11^(0, 0) = g 2 



{/i 2 - m 2 ) 1/2 /i 



^ In 



p + (p 2 



7T 7T 

The ultimate expression for the electric mass has the form 



m 



2\l/2 



m 



m 



el 



g 2 T 2 



N r . 



+ En / (o,o) 



/=i 



(20) 
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In this approximation, the effect of the instanton liquid is completely accounted for by the quark 

mass dynamically generated in the instanton medium. With such definition of mass, the formula (|TU|) 

1 2 
at p = and T^O should be modified. The coefficient ^ at Nf should be replaced by — y. However, 

this replacement has only a little effect; self-consistency of our calculations will be discussed below. 

Using the integral f|T7|) . which is also valid for the caloron solution, we derive the expression for 
the distribution of pseudoparticles: 



d(p; /i, T) = d(p; 0, 0) e 



-r?{p,T) p 2 



2 7T Z 



rp'2 



Nr. 



E n/ (o,o) 



(21) 



The one-loop quark contribution to the instanton action at zero temperature, finite chemical poten- 
tial, and w ^ was studied in detail in [H] (see also [T5] . [T5]). These studies make it possible to 
improve our description, however, we work within the approximation (I2T1) and, moreover, we consider 
the limit of massless quarks. A self-consistent calculation for the quark with dynamically generated 
mass can be the subject of a separable study. 

Necessary modifications in the variational principle are as follows. It was revealed that only 
the distribution function d(p; p, T) of pseudoparticles changes, whereas the repulsion interaction 
Uint between pseudoparticles remains as before. Similar to the case of instantons, we introduce the 
parameter v satisfying the relation 

(22) 



P 1 



instead of (jHJ). Since the instanton liquid density is greater than zero, a new limitation on the average 

size of pseudoparticle emerges pA < • If this limit is smaller than the limit descussed above, 
then it must be the starting point for the determination of the equilibrium size of pseudoparticles by 
the bisection method. The derivative of the function (3 with respect to the density of the instanton 
liquid can be determined from the relation f[2"2"j) . The result is 



n d/3 
(3 dn 



v—r] z p z 



4/3 



(23) 



it should be substituted in Eq. ( TTOl that determines the saddle-point. The integral ( |T9l) should be 
evaluated numerically because it cannot be calculated analytically at arbitrary temperatures even 
though the quark mass equals zero. Thus we are ready to determine the parameters of the instanton 
liquid everywhere over the p - T plane. For simplicity, the calculations are performed at zero quark 
masses. We neglect the light-quark contribution to the respective determinants [T7j (see also [18]). 
We also disregard a possible temperature molecular behavior of instanton-anti- instanton pairs |19j . 

The results of the calculations are shown in Fig. 1 by the lines of constant density. The instanton 
liquid density is plotted in Fig 2. versus the temperature (at zero chemical potential) and versus the 
chemical potential (at zero temperature). Though the conventional natation for the instanton liquid 
density at nonzero temperature is n = TN/V3, we use the label n which is more simple. At T 7^ 0, 
and p = 0, our results coincide with the results obtained in [UJ and [9]. It sould be noted that our 
results are consistent with recent calculations on a lattice at finite temperatures [20], [21], where a 
rapid decrease of the chromoelectric components in the respective correlation functions was found. 
In our model, such suppression is due to the term rj m ez ^4^4 m the effective action; with neglect of 
this factor, the chromoelectric and chromomagnetic correlators coincide. From this point of view, our 
calculations may seem inconsistent. We use the caloron solution (1131) . which is symmetric under an 
interchange of chromoelectric and chromomagnetic fields. However, the caloron components manifest 
themthelves in the observables differently because of the anisotropy of the weight function. In fact, 
our method of taking the gluon mass term into account is consistent only in perturbation theory. In 
a complete study, one must find an analogue of the solution ([TBI for the effective Lagrangian with the 
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Figure 1: Lines of equal density of the instanton liquid in the temperature-chemical potential plane. 
Curve 1 corresponds to the density n = 0.75 n , where n is the density at zero temperature and 
chemical potential. Also shown are the densities from (curve 2) n — 0.5 tiq to (curve 6) n — 0.1 n®. 
Curves 3-5 correspond to intermediate densities at intervals of 0.1. 

0.8 i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 




T/A, M /A 

Figure 2: Instanton liquid density versus (curve 1) temperature and (curve 2) chemical potential. 

gluon mass generated for the chromoelectric field and gain a self-consistent description of ensemble 
of pseudoparticles in the long- wave approximation [3] . 

It is of interest that the data on correlation functions for cooled configurations [20] are fitted 
well by the instanton ensemble [22] • In so doing, the contribution of the terms of the second order 
in the instanton liquid density (~ n 2 ) is in excellent agreement with the effect of the standard 
instanton ensemble with the respective admixture of the perturbative component everywhere over the 
distance range chosen for a fit [23]. This agreement indicates that the confining component is absent 
from the lattice configurations isolated by cooling. It is surprising because lattice simulations with 
cooling were aimed at the searches for a long- wave confining component. However, an interpretation 
of lattice simulations at finite temperature presents difficulties because it is not clear what scale 
corresponds to the configarations used for the measurements. The magnitude of deformation of the 
chromoelectric component of the solution for the effective Lagrangian with the mass term is also 
poorly known. The scale of lattice configurations can, in principle, be estimated using the scale at 
which the chromoelectric field decrease since only this scale has emerged in our calculations. 

In conclusion we note that, though we used only a rough approximation, the most important 
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features of the behavior of the instanton liquid density (gluon condensate) in the medium have been 
revealed. The lines of equal density are markedly extended along the \i axis because, according to 
the formula (1201) . the most substantial gluon component of screening vanishes at small temperatures. 
Typical values of T and fi at which the effects of the medium become significant are related to each 

other by the formula ^ c ^i^J-l^l (T/A) 2 ~ Nf ^ 1 ; which leads to a plausible coefficient 

of oblongness along the \i axis 



(at N c = 3 and Nf = 2). A fall in density evaluated with allowance for the dynamically generated 
quark mass should begin at a greater value and be more steep. The reason is that, at chemical 
potentials less than the quark mass, the quark contribution to screening is reduced. This gives rise 
to formation of a plateau and concentration of the lines of equal density. The dependence of the 
dynamical quark mass on the momentum u is significant at small temperatures leading to a decrease 
of screening approximately by a factor of two [13] . 

We are grateful to A.E. Dorokhov for helful discussions. 
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